Ab initio many-body calculations of deuteron- 4 He scattering and 6 Li states 
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We extend the ab initio no-core shell model/resonating- group method (NCSM/RGM) to 
projectile-target binary-cluster states where the projectile is a deuteron. We discuss the formal- 
ism in detail and give algebraic expressions for the integration kernels. Using a soft similarity- 
renormalization-group evolved chiral nucleon-nucleon potential, we calculate deuteron- 4 He scatter- 
ing and investigate 6 Li bound and unbound states. Virtual three-body breakup effects are obtained 
in an approximated way by including excited pseudo-states of the deuteron in the calculation. We 
compare our results to experiment and to a standard NCSM calculation for 6 Li. 
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I. INTRODUCTION 

Ab initio many-body calculations of nuclear scattering 
and reactions pose a challenge to nuclear theory. For A — 
3 and 4 nucleon systems, the Faddeev [l| and Faddeev- 
Yakubovsky Q as well as the hyperspherical harmonics 
(HH) or the Alt, Grassberger and Sandhas (AGS) 
[J] methods are applicable and successful. For systems 
with more than four nucleons, only very few approaches 
presently exist among which the Green's function Monte 
Carlo method applied recently to the calculation of the 
n- 4 He scattering Q. 

Recently we combined the resonating-group method 
(RGM) and the ab initio no-core shell model 

(NCSM) [l2j. into a new many-body approach fl3l - fll| 
(ab initio NCSM/RGM) capable of treating bound and 
scattering states of light nuclei in a unified formalism, 
starting from the fundamental inter-nucleon interactions. 
So far, applications have been limited to the description 
of projectile-target scattering where the projectile is a 
single nucleon. In particular, we first studied the n- 3 H, 
n- 4 He, p- 3 ' 4 He, n- 10 Be, scattering processes [H,[II[ and 
later also the n- 7 Li, p- 7 Be scattering as well as nucleon 
scattering on 12 C and 16 O 15|. In the present paper we 



extend the formalism to the case of a two-nucleon projec- 
tile and perform calculations for deuteron- He (or d-a) 
scattering. Simultaneously, we investigate the d-a bound 
state and compare our results to a standard NCSM calcu- 
lation for 6 Li. It should be emphasized that the present 
formalism is general and applicable to any other target 
nucleus, i.e. to any deuteron-nucleus system. 

The deuteron is weakly bound and can be easily de- 
formed. Its polarization and virtual breakup cannot be 
neglected even at very low energies. A proper treat- 
ment of these effects requires the inclusion of three- 
body continuum states: neutron-proton-nucleus. This 
is very challenging. Even though the extension of the 
RGM formalism to include three-body clusters is feasi- 
ble [IH, [l7| , in this first application we limit ourselves to 
two-body clusters only and approximate virtual three- 
body breakup effects by discretizing the continuum with 
excited deuteron pseudo-states. 

Deutcron- 4 He scattering was investigated within the 



binary-cluster RGM formalism in the past jl8H2C| . How- 
ever, the present investigation is the first that uses accu- 
rate nucleon-nucleon (NN) interactions (i.e. such that fit 
the NN phase shifts with high precision) and many-body 
cluster wave functions obtained consistently from the 
same Hamiltonian. We do not fit or adjust any parame- 
ters, rather we systematically investigate the convergence 
of our results with respect to the size of the harmonic os- 
cillator (HO) basis used to expand the cluster wave func- 
tions and localized parts of the RGM integration ker- 
nels as well as with respect to the number of deuteron 
pseudo-states and/or 4 He excited states included in the 
calculation. We compare our results to a standard ab 
initio NCSM calculation for 6 Li that uses the same NN 
potential. In this study, we employ a similarity rcnormal- 
ization group (SRG) (2l|, [H] evolved chiral N 3 LO JVJV 
potential [23| (SRG-N 3 LO) that is soft enough for us to 
reach convergence within about 12— 14M7 HO excitations 
in the basis expansion. 

In Sect. [H] we briefly overview the general features 
of the NCSM/RGM formalism and present for the first 
time algebraic expressions for the NCSM/RGM integra- 
tion kernels when the projectile nucleus has mass number 
a = 2. The matrix elements of the norm kernel are given 
in this section, while those of the Hamiltonian kernel are 
presented in Appendix [AJ In Sect. IIII1 we discuss our 
results for d-a scattering and bound-state calculations. 
We show the calculated phase shifts and cross sections 
and compare the deuteron- 4 He results to 6 Li ab initio 
NCSM calculations with the same Hamiltonian. Conclu- 
sions and outlook are given in Sect. [TV] 



II. FORMALISM 



In the present paper we a pply the NCSM/RGM for 



malism introduced in Rcf. [14} to the description of 
deuteron-nucleus collisions. While the derivation of the 
integration kernels was specialized for projectile-target 
basis states with a single- nucleon pro j ectile , the theoret- 
ical framework presented in Ref. [lj] is general and fully 
applicable to the present case. In this section we briefly 
revisit the NCSM/RGM formalism and provide algebraic 
expressions of the integration kernels for the specific case 
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of a two- nudcon projectile. 

Following the notation of Ref. [l4| , the wave function 
for a scattering process involving a two-nucleon projectile 
and a target nucleus can be cast in the form 



are the Hamiltonian and norm kernels, respectively. Here 
E is the total energy in the center-of-mass (cm.) frame, 
and H is the intrinsic A-nucleon microscopic Hamilto- 
nian, for which it is useful to use the decomposition, e.g.: 



(1) 



through an expansion over binary-cluster channel-states 
of channel spin s, relative angular momentum £, total 
angular momentum J, parity 7r, and isospin T, 
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The above basis states are uniquely identified by the 
channel index v = {A-2 a^Tx; 2a 2 IpT 2 ; si}. The 
internal wave functions of the colliding nuclei contain 
A— 2 and 2 nucleons (A>2), respectively, are antisym- 
metric under exchange of internal nucleons, and depend 
on translationally invariant internal coordinates. They 
are eigenstates of H^a-2) an d H( 2 ), the (A— 2)- and two- 
nucleon intrinsic Hamiltonians (li, 7Tj, Ti and a>i denote 
respectively spin, parity, isospin and additional quantum 
numbers of the i-th cluster) . The clusters centers of mass 
are separated by the relative vector {fi being the position 
vector of the i-th nucleon) 



H = T rel (r) + V re] + V c ir) + H (A _ 2) + H, 



(8) 



Further, T rc i(r) is the relative kinetic energy and V re i 
is the sum of all interactions between nucleons belong- 
ing to different clusters after subtraction of the average 
Coulomb interaction between them, explicitly singled out 
in the term Vc{r) = Z\ v Z 2v e 2 / f r ', where Z\ v and Z 2v are 
the charge numbers of the clusters in channel v. 
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(9) 



The Vij interaction consists of the strong and Coulomb 
part. Thanks to the subtraction of Vc(r), the overall 
Coulomb contribution presents a r~ 2 behavior, as the 
distance r between the two clusters increases. Therefore, 
V ro i is localized also in presence of the Coulomb force. 
In this paper, we limit our calculations to the use of a 
two-nucleon interaction only, but the formalism can be 
generalized to include the three-nucleon interaction in a 
straightforward way. 
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In Eq. (|2|) , the residual anti-symmctrization for exchange 
of nucleons pertaining to different clusters is guaranteed 
by the anti-symmetrizer for the {A— 2,2) partition 



A v = A(A- 2 ,2) 

A- 2 A A-2 

1 - E E p * + E Ka-iPj.a 
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,(4) 



A. Norm kernel 

For definitions and details regarding the derivations 
outlined in this and the next section we refer the inter- 
ested reader to Sees. II.C.l and II.C.2 of Ref . [3. 

Because the wave functions of both (A-2)-nucleon and 
two-nucleon clusters arc anti-symmetric under exchange 
of internal nucleons, the norm kernel ([7]) for the same, 
(A— 2, 2), mass partition in both the initial and final state 
can be written as 



where C is the normalization constant 



The 



A(A-l) 

unknown relative-motion wave functions g„ T (r) can be 
determined by solving the many-body Schrodinger equa- 
tion in the Hilbert space spanned by the basis states 
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where 



dr r 



U j ; u T {r',r)-ENZ: v T {r',r) 



9 J J T {r) 



n J JJ(r',r) 



Ay' HA 



■ IT 



-- 0, 

(5) 

(G) 
(7) 



(a) 



(b) 




FIG. 1: (Color online) Diagrammatic representation of: (a) 
"direct", (b) "one-nucleon-exchange" and (c) "two-micleon- 
exchange" components of the norm kernel. The groups of 
circled lines represent the (A— 2)- and two-nucleon clusters. 
Bottom and upper part of the diagram represent initial and 
final states, respectively. 
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where |$^ n T ) is the translationally-invariant HO channel 
state introduced in Eq. (22) of Ref. [l4j . here for clarity: 



If 



(|A-2 ai / 1 7ri T 1 )|2a 2 / 2 ^T 2 ) 



x Y<? (fA-2,2) 



R n i{r A - 2a ) . (12) 



Three terms contribute to the norm kernel (JTOJ): A di- 
rect term, in which initial and final states are identical 
(corresponding to diagram (a) of Fig. [1]); a one-nucleon 
exchange term, corresponding to diagram (b) of Fig. [T] 
and, finally, a two-nucleon exchange term, corresponding 
to diagram (c) of Fig. [TJ 

In this paper, the localized parts of the integration 
kernels ([6]) and (O are obtained in two steps. First, ma- 
trix elements of translationally-invariant operators (for 

the norm P A -i,A and Pa-2,aPa-3,a-i) are calculated 
in the Slater-determinant (SD) basis, in which the eigen- 
states of the (A— 2)-nucleon fragment are expanded in 
HO Slater determinants: 



1$ 



/SD 



|A-2a 1 / 1 T 1 ) SD |2a 2 / 2 T 2 

-1 (J*T) 



R,u(RgL). (13) 



Second, the corresponding translationally-invariant ma- 
trix elements on the basis (jT2"j) are recovered throu gh a 
transformation as described in Sect. II. C. 2 of Ref. [14j . 
Eq. (32). Here we remind that the eigenstates of the 
(A— 2)-nucleon fragment in the SD basis, 

(fv TA-iov ■ ■ aA-2Tr • • T J 4- 2 |^-2ai/ 1 ri Ti)sD, (14) 

are related to the transaltionally-invariant eigenstates by 
the expression 

|A-2a 1 / 1 T 1 ) SD = \A-2 ai hT x ) ^(R^) , (15) 

and the cm. coordinates introduced in Eqs. (|13j) and 



15|) are given by: 



■A-2 



=A-1 



(16) 




The calculation of matrix elements in the basis ([13)) is 
most efficiently achieved by first performing a transfor- 
mation to a new SD basis: 



x(-i) h+J+e+e ^sii 2 s 2 jaj b Ll b 

x (nJ a n b £ b L ab \n£n 2 £ 2 L ab ) d=1 

x (n 2 £ 2 s 2 I 2 T 2 \2a 2 I 2 T 2 ) \$ J k J)sd (17) 

where the sum runs over the quantum numbers n 2 , £ 2 , s 2 , 
n a ,i a ,ja, nb,£b,jb, L ab , and /, (n 2 £ 2 s 2 I 2 T 2 \2 a 2 I 2 T 2 ) 
is the projectile wave function expanded in the 
relative-coordinate HO basis, s = s/2s+l etc., and 
(nal a nbhLab\n£n 2 £ 2 Lab}d=i indicates an HO bracket 
for two particles with identical masses. In addition, 
we introduced the cumulative quantum number n ab = 
{A-2 a\hT\; n a £ a j a \;n b £ b j b \;IT 2 } and the new SD 
channel states: 



K7>sd 



(IT2) 
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(18) 



Using the basis states of Eq. (fTBf to evaluate the ma- 
trix elements of the transposition operators appearing in 
Eq. (fT0|) results in the following expressions: 
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(20) 



where the indexes a and 6 represent the sets of single- 
particle quantum numbers {n a (. a ja\} and {nbibjb^}-, re- 
spectively, such that at = a t . i etc., a' and 6' are 

analogous indexes associated with the primed quantum 
numbers, n' ab = {A-2 a[I^T{; a ';b '; / T 2 '} , and, finally, 
Vjmim, = (-ir m+ 5-""a tt(j _ ra i_ rat . In addition, we 
note that Eqs. (fT9|) and (|20|) depend on the one- and two- 
body density matrix elements (OBDME and TBDME), 



respectively, of the target nucleus. 



B. Hamiltonian Kernel 

The Hamiltonian kernel © for the same, (A— 2,2), 
mass partition in both the initial and final state can be 
cast in the form 



v ■ r 



A 
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where the potential kernel is defined by 
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Clearly, for the (A— 2,2) partition the potential kernel norm kernel. We identify five separate terms corrcspond- 
presents a much more complicated expression than the 



5 




(a) (b) (c) 




(d) (e) (f) (g) 

FIG. 2: (Color online) Diagrammatic representation of: "di- 
rect potential", (a)-(b), and "exchange-potential", (c)-(g), 
components of the potential kernel. The groups of circled 
lines represent the (A— 2)- and two-nucleon clusters. Bot- 
tom and upper part of the diagram represent initial and final 
states, respectively. 

ing to the nine diagrams presented in Fig. [2] The first 
"direct-potential" term on the right-hand side (rhs) of 



Eq. (f23|) corresponds to diagrams (a) and (b) , the second 
term corresponds to diagram (c) , while diagrams (d) and 
(e) represent the third term. The last two terms are then 
depicted schematically by diagrams (f) and (g), respec- 
tively. Matrix elements of each of these terms in the basis 
(fT8|) are given in Appendix El The first two terms, (|Alj) 
and (|A2[) . depend on the OBDME of the target nucleus; 
the second two terms, (|A3|) and (|A4[) . depend on the TB- 
DME of the target nucleus; and, finally, the last term, 
(|A5|) . depends on the three-body density of the target 
nucleus. The three-body density matrix elements can be 
obviously re-coupled in different ways. Here, we selected 
a particular angular-momentum coupling that results in 
the simplest expression for our Hamiltonian kernel matrix 
element. In general, it is a challenge to compute three- 
body density matrix elements, in particular due to their 
rapidly increasing number in the multi-major shell basis 
spaces. However, since in the present paper we focus on 
the A = 6, G?- 4 He and 6 Li systems, we can take advantage 
of the completeness of the (A— 5)-body eigenstates and 
re-write the expression (|A5[) in the form 
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(24) 



where the sum runs over the quantum numbers /3, /g, Tp 1 
d = n d i d j d \ etc., e,e', J de ,T de , X,Y,t x , and t y . For 
the present case of A = 6, the states \A— 5 (3 IpT^sn 
reduce to the HO single particle states In^/^j^i) and 
the reduced matrix elements in Eq. (f2"4"]l that involve 4 He 
eigenstates are straightforward to calculate. 

We also note that the terms (|A1[) and (|A3|) are 
symmetric, while the remaining ones, (| A2|) . (|A4[) and 
(|A5p . are not. Therefore, we introduce a Hermitizcd 
NCSM/RGM Hamiltonian, as discussed in detail in 
Ref. [TI, using AHA = \{A 2 H + HA 2 ) (see Eq. (42) 
in RefTpj). 



III. APPLICATION TO THE DEUTERON- 4 HE 
SYSTEM 



The deuteron-nucleus formalism presented in the pre- 
vious section is completely general. The simplest system 
to which it can be applied is deuteron- 4 He for two rea- 
sons. First, the complicated calculation of the target 
three-body density needed to compute the last term on 
the right-hand side of Eq. ([^l . given by Eq. (|A5[) . be- 
comes straightforward for for A = 5 and 6 (that is for 
3 H, 3 He and 4 He targets) using the completeness of the 
(A— 5)-nucleon eigenstates as demonstrated in Eq. (|24p . 
Second, the 4 He nucleus is tightly bound with its first 
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Eg. e . [MeV] 


2 H 


4 He 


6 Li (NCSM/RGM) 


6 Li (NCSM) 


Calc. 


-2.20 


-28.22 


-32.25 


-32.87 


Expt. 


-2.22 


-28.30 


-31.99 


-31.99 



TABLE I: Calculated g.s. energies of 2 H, 4 He, and 6 Li ob- 
tained by using the SRG-N 3 LO NN potential with A = 1.5 
fm -1 are compared to the corresponding experimental val- 
ues. The NCSM calculations for 2 H, 4 He and 6 Li were per- 
formed in -/V max = 12, 12 and 10 basis space, respectively. The 
NCSM/RGM calculation included 4 He and 2 H ground states 
and 7 deuteron pseudo-states in each of the 3 Si- 3 Di and 3 D3- 
3 G3 channels, as well as 5 pseudo-states in the 3 Z>2 channel. 
The HO frequency of hQ = 14 MeV was used. 



excited state at E x ss 20 MeV. A solution of the coupled- 
channel equations ([5]) obtained limiting the target-states 
to the ground state (g.s.) is already a very good approx- 
imation for the d-a system. 

To test the formalism, we use a soft SRG-N 3 LO NN 
potential with evolution parameter A = 1.5 fin -1 . With 
this low value of A, our calculations reach convergence at 
A max ps 12. We note, however, that a somewhat higher- 
value of A would result in a better agreement with exper- 
imental data as discussed later. We benchmark our d-a 
NCSM/RGM results with standard NCSM calculations 
for 6 Li. Any differences can then be attributed to miss- 
ing degrees of freedom rather than to the model space 
truncation. 

Our calculation starts with the NCSM diagonalization 
of the Hamiltonian in the N max hfl HO basis for d and 
4 Hc. Obtained eigenenergies and eigenfunctions serve 
then as input in Eq. (JSJ) - First, one-, two-, and three- 
body densities are calculated from the 4 He wave func- 
tions, then the integration kernels are calculated. The 
localized parts of the integration kernels are expanded in 
the same A^ max (or N max + 1 depending on parity) HO 
basis space as the cluster cigenstates. The same HO fre- 
quency is used in all calculations. The wave functions 
of the d-a relative motion are found by solving (JS|) with 
either bound-state or scattering-state boundary condi- 
tions by means of the microscopic R-matrix method on 
a Lagrange mesh [24| (with additional details given in 
Ref. |l3, Sect. II. F). 



A. Bound-state calculations 

Our results for the ground states of deuteron, 4 He and 
6 Li are presented in Table |TJ The convergence of the 
NCSM calculations can be judged from Fig. [3l where we 
show both absolute and excitation energies of 6 Li as well 
as the d+a threshold. The 4 He convergence is excellent 
and that of 6 Li very good. The 6 Li excited states are 
resonances, but within the NCSM calculation they are 
approximated by eigenstates expanded in the HO basis. 
The 4 He is slightly under-bound while the 6 Li is over- 
bound by about 0.9 MeV due to the choice of a low A 
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FIG. 3: (Color online) Absolute (top panel) and excitation 
(bottom panel) energies of 6 Li calculated within the NCSM 
compared to experiment. The dashed lines in the top panel 
indicate the calculated NCSM (A max = 8, 10, 12) and exper- 
imental d-a thresholds. The SRG-N 3 LO NN potential with 
A = 1.5 fm -1 and the HO frequency of hfl = 14 MeV were 
used. 



value and the neglect of the SRG-induced three-body in- 
teraction]^ . We note that in the NN-only calculations 
of Ref. [25|], selecting A«2 fm" 1 results in 4 He and 6 Li 
binding energies closer to experiment. The excited states 
are correctly ordered except for the reversal of the 2+1 
and the 1^0 states. The splitting of the 2 + and the 
3 + states is underestimated, a sign of weak spin-orbit 
interaction, most likely due to the neglect of the initial 
chiral three-nucleon interaction. One more feature to no- 
tice is the drop of the 2 + and 1^0 excitation energies 
with increasing N max . This is a consequence of the fact 
that these states are broader resonances compared to the 
3+0 or the T = 1 states. 

The variation with respect to N max of the 6 Li g.s. 
energy calculated within the NCSM/RGM is similar to 
that of the 4 He. The N max — 12 model-space is suf- 
ficient to reach convergence. However, the deuteron is 
weakly bound: Its polarization and virtual breakup can- 
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FIG. 4: (Color online) Dependence of the 6 Li g.s. energy on 
the number of deuteron pseudo-states Ad* in the 3 <Si- 3 -Di, 
3 Z?2 and 3 D3- 3 Gs channels included in the calculation. Di- 
amond symbols: All three channels included. Triangle sym- 
bols: Only 3 Si- 3 L>i channel included. The NCSM results 
(solid line) and the 2 H+ 4 He threshold energy (dashed line) 
are also shown. Details of the calculations are described in 
the caption of Table |U 



not be neglected. This is demonstrated in Fig. |4] The 
NCSM/RGM calculation limited to the deuteron ground 
state binds 6 Li by only about 200 kcV contrary to the 
NCSM result of 2.4 MeV using the same Hamiltonian. To 
include deuteron polarization and virtual breakup prop- 
erly one would have to extend the NCSM/RGM formal- 
ism to a three-cluster basis: n-p-a. This is quite chal- 
lenging. In the present work, we discretize the contin- 
uum by including excited deuteron pseudo-states in the 
NCSM/RGM coupled channel equations. The pseudo- 
states are obtained in the NCSM diagonalization. In Ta- 
ble HH we present the pseudo-state energies obtained in 
the iV max = 12 basis. The 6 Li g.s. convergence with re- 
spect to the number of d pseudo-states Ad* included in 
the calculation is shown in Fig. 01 The S'-wave dominated 
(odd Ad*) pseudo-states in the 3 Si- 3 D\ channel have a 
quite dramatic influence on the 6 Li binding energy. The 
pseudo-states in the 3 i?2 and 3 D2,- 3 G^ channels are less 
important for the ground state but have a significant ef- 
fect on the 2 + and the 3 + 6 Li resonances. By including 7 
or 9 pseudo-states, we reach convergence with respect to 
the number of d* in the channels considered here (for the 
g.s. energy, the Ad* = 7 and 9 results are within 30 keV 
of eachother) . Still, the 6 Li NCSM calculation contains 
more correlations as it produces a lower g.s. energy by 
about 600 keV or 2%. This can be seen from Table U 
and Fig.|H The missing correlations in the NCSM/RGM 
calculation most likely include excitations of the 4 He (of 
which we have included here only the ground state), as 
well as deuteron excitations in other channels. To es- 
timate to which degree excited states of the 4 He target 
would influence our NCSM/RGM results for the 6 Li g.s. 
energy, we performed a calculation that included g.s. and 
first excited + state of 4 He as well as deuteron g.s. and 



F [MoVl 


3 Si- 3 Di 


3 Do 


3 D?- 3 G? 


g.s. 


-2.20 






1 * 






7 fil 


2'* 


7 RQ 


1 8 81 
10.01 


1 ^ 79 


Q* 
•i 


15.20 


35.05 


19.32 


4* 


19.74 


57.28 


33.13 


5* 


31.90 


87.88 


36.47 


6* 


37.60 




57.01 


T 


55.95 




60.13 



TABLE II: Calculated H g.s. and pseudo-state energies ob- 
tained using the SRG-N 3 LO AW potential with A = 1.5 
fm - in the A max = 12 basis space and HO frequency of 
= 14 MeV. In our largest calculations, 7 pseudo-states 
were included in the two coupled channels and 5 pseudo-states 
were included in the 3 Di channel. 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




[fm] 



FIG. 5: (Color online) Ground-state wave function of 6 Li 
as a function of the separation between deuteron and 4 He 
clusters. Dominant 5-wave (solid lines) and D-wave (dashed 
lines) components are shown. The symbol d* (d'"*') denotes 
the first (n-th) deuteron pseudo-state of the 3 Si- 3 Di chan- 
nel, while d"" denotes the lowest deuteron pseudo-state of the 
3 Dz- 3 G3 channel. Details of the calculation are described in 
the caption of Fig. U] 

7 d* pseudo states in the 3 S\- 3 Di channel. The binding- 
energy increased by 71 keV compared to the calculation 
with the same number of deuteron states but only the 
ground state of 4 He. This is a no n- negligible effect and, 
based on our previous study of the n- 4 He system [l4|; it 
is plausible that the addition of the next five or so lowest 
excited states of 4 He would provide an extra ~ 500 kcV 
of binding. Unfortunately, such a calculation is currently 
out of reach. 

Finally, we note that contrary to the NCSM calcu- 
lations, the NCSM/RGM bound state has the proper 
asymptotic behavior of a Whittaker function with re- 
spect to the d+a threshold. In Fig. [SJ we can see that 
the S'-wave extends well beyond 10 fm. This plot of the 
NCSM/RGM wave function can be compared with Fig. 
6 of Ref. [2(|. There the overlap functions [g(r)/r] of 
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the 6 Li ground state with <i+ 4 He cluster states obtained 
within the standard NCSM vanish beyond about 8 fm. 



B. Scattering calculations 

By solving the NCSM/RGM coupled-channel equa- 
tions ([5]) for positive energies, we obtain the wave func- 
tions of the relative motion of the clusters and the scat- 
tering matrix for each considered JT channel. The scat- 
tering matrix can then be used to calculate cross sections 
and other observables. 

In Figs. |B]|n]we present our calculated diagonal S- and 
D-wave phase shifts. First, we study the phase-shift con- 
vergence with respect to the size of the HO basis expan- 
sion for the cluster wave functions and localized parts 
of the integration kernels. In Fig. HI we show phase 
shifts obtained respectively in the N max = 12 (solid line), 
iV ma x = 10 (dashed line), and iV max = 8 (dotted line) 
model spaces. The curves in the top panel include only 
the ground states of d and 4 He, while the middle and 
bottom panels show results including up to 7 deuteron 
pseudo-states in the 3 <Si- 3 -Di and 3 D 2 and 3 -D 3 - 3 G 3 chan- 
nels. The N max = 10 and N max = 12 lines are on top of 
each other in the J^T = 1+0 channels while some small 
change in the phase shifts is still visible in the 2+0 and 
the 3 + channels. These differences become smaller in 
calculations with the pseudo-states. Overall, the conver- 
gence is satisfactory. At this stage an -/V max = 14 calcula- 
tion would be computationally very challenging. Figure[7] 
demonstrates the phase shift convergence with respect 
to the number of pseudo-states included in the coupled- 
channel NCSM/RGM equations. It is clear that, similar 
to the bound-state calculation, for the d* channels con- 
sidered here convergence is reached with 7 pseudo-states. 
The relative contribution of pseudo-states from the three 
d* channels considered here can be judged from Fig [3] 
The 3 Si- 3 Di pseudo-states affect all S and D waves. 
On the contrary, the 3 Z?2 and the 3 Z?3- 3 G3 pseudo-states 
have a considerable effect only on the 3 £>2 and 3 D$ waves. 
Note that the solid lines in Fig. [8] correspond to the dot- 
ted lines in Fig. [7] 

Our calculated diagonal S- and D-wave phase shifts 
are compared to the phase shifts extracted from experi- 
mental data in Refs. [27j and [28[ in Fig. [§] The calcula- 
tion corresponds to the largest basis space (N amx = 12) 
and the highest number of deuteron pseudo-states that 
we employed in this work. Our S-w&ve and 3 Z?3-wave re- 
sults compare well with the experimental data. However, 
the 3 Di and in particular the 3 D 2 phase shifts overesti- 
mate the experimental ones. The position of our calcu- 
lated 2 + resonance is below the experimental one by 
almost 1 MeV. The splitting between the D-waves is 
underestimated. Clearly, the strength of the spin-orbit 
interaction in the calculation is smaller than it should 
be. This is most likely due to the neglect of the three- 
nucleon forces in our calculations, those induced by the 
SRG transformation and, more importantly, the initial 
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FIG. 6: (Color online) Calculated d- 4 He S- and D -wave phase 
shifts. Solid, dashed and dotted lines correspond to the 
■/Vmax = 12, 10 and 8 basis sizes, respectively. Results in the 
top panel were obtained considering only the ground state of 
the deuteron projectile. In the middle panel, calculations in- 
corporate 7 additional deuteron pseudo-states in the 3 Si- 3 Di 
channel. In the bottom panel, up to 7 deuteron pseudo-states 
were included also in the 3 D2 and 3 Z?3- 3 G3 channels. The 
SRG-N 3 LO NN potential with A = 1.5 fm" 1 and the HO 
frequency of = 14 MeV were used. 
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FIG. 7: (Color online) Calculated d- 4 He S- and D-wave 
phase shifts. Dependence on the number of included deuteron 
pseudo-states in the 3 Si- 3 Z?i, 3 Di and 3 Ds- 3 G3 channels: 
Solid, dashed, dotted, and thin dotted lines correspond to 
up to 7, 5, 3, and 1 pseudo-states in each channel, respec- 
tively. The SRG-N 3 LO NN potential with A = 1.5 far 1 , 
the iV m ax = 12 basis size and the HO frequency of HQ = 14 
MeV were used. 




FIG. 8: (Color online) Calculated d- 4 He S- and D -wave 
phase shifts. Influence of deuteron pseudo-states from differ- 
ent channels. Solid lines correspond to calculations with one 
deuteron pseudo-state in each of the 3 Si- 3 Di, 3 D 2 and 3 D3- 
3 G3 channels. Dashed lines identify results obtained with only 
one deuteron pseudo-state in the 3 Si- i Di channel. Dotted- 
line curves are the solution of calculations including only the 
ground state of the deuteron. The SRG-N 3 LO NN potential 
with A = 1.5 fm -1 , the iV max = 12 basis size and the HO 
frequency of hQ = 14 MeV were used. 



chiral EFT three-nucleon interaction. Our calculated P- 
and F-wave phase shifts are presented in Fig. [TUJ While 
the _F-waves monotonically increase, the P-waves exhibit 
more structure and, in particular, the 3 Pq changes sign 
and becomes negative beyond the center-of-mass energy 
E kin > 2 MeV. 

In Fig. [TTJ we compare our calculated differential cross 
section to the experimental data of Rcfs. (2|| and [3(| 



FIG. 9: (Color online) Calculated d- 4 He S- and D -wave phase 
shifts compared to experimental data from Refs. [2tJ (circles) 
and [111 (diamonds). Up to 7 deuteron pseudo-states were 
included in each of the 3 S l i- 3 Di and 3 Ds- 3 G3 channels and 
5 pseudo-states in the 3 D 2 channel. The SRG-N 3 LO AW 
potential with A = 1.5 far 1 , iV m „ = 12 basis size and HO 
frequency of Ml = 14 MeV were used. 




1 2 3 4 5 6 7 



FIG. 10: (Color online) Calculated d- 4 He P- and F-wave 
phase shifts. Details of the calculation are given in the caption 
of Fig. [9] 



for four deuteron laboratory energies in the range Ed rj 
3 — 12 MeV. Our calculation over-predicts the measured 
cross section at E c i = 2.94 MeV, most likely a conse- 
quence of the incorrect position of the calculated 2+0 
resonance, see Fig[5J However, for the intermediate en- 
ergies, Ed = 6.97 and 8.97 MeV, the agreement with 
the measured data is reasonable. At Ed = 12 MeV the 
differences become larger. We also note that our calcu- 
lated cross section underestimates the data in the range 
of # c .m. ~ 20 — 45 deg. It should be kept in mind that 
in our calculations the deuteron breakup is accounted for 
only by using the pseudo-states rather than as a three- 
body final state. To shed light on the influence of the 
pseudo-states on the cross sections, we show in Fig. Q2] 
results for Ed = 6.97 MeV obtained with pseudo-states 
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FIG. 12: (Color online) Influence of deuteron pseudo-states 
from different channels on the d- 4 He differential cross section 
at the deuteron laboratory energy of 6.965 MeV. The experi- 
mental data (symbols) are from Ref. [2jj. The solid line cor- 
responds to the calculation described in Fig. [9] The dashed- 
dotted line indicates the results obtained with 7 deuteron 
pseudo-states in the 3 Si- 3 Di and 3 Z?3- 3 G3 channels. The 
dashed-line curve is the solution of the calculation with 7 



deuteron pseudo-states only in the 
waves up to J=6 were included. 



'Si- Di channel. Partial 



FIG. 11: (Color online) The d- 4 He differential cross section 
at the deuteron laboratory energies of 2.935, 6.965, 8.971 and 
12 MeV. The experimental data (symbols) are from Ref. [2Sj 
and [13]. The calculations (lines) are as described in Fig. [9] 
Partial waves up to J— 6 were included. 



in different channels. With the pseudo-states only in the 
3 Si— 3 Z?i channel, the cross section is not well described 
beyond 50 deg. The inclusion of the pseudo-states in 
the 3 D3- 3 G3 channel improves the agreement with the 
data somewhat. By adding the pseudo-states in the 3 D 2 
channel, the agreement with the data beyond 50 deg is 
quite reasonable but at the forward angles, from 20 deg 
to 45 deg, the agreement with the data is spoiled. The 
NN interaction that we employed is not the optimal one 
as explained earlier: A low value of A = 1.5 fm -1 was 
selected to facilitate a fast convergence and a straight- 
forward comparison to the standard NCSM calculation. 
As seen in Fig. |H1 this potential overestimates the 2 + 
D-wave phase shifts compared to the data. The pseudo- 
states from the 3 Z?2 channel enhance this overestimation 
as visible in Fig. [5] This is the likely cause of the wors- 
ening of the cross-section agreement with the data at 
forward angles when the Z D 2 pseudo-states are added to 
the NCSM/RGM basis. 

The g.s. and the T = resonance-state energies ob- 
tained within the NCSM and the NCSM/RGM are com- 
pared to each other and to experimental values in Fig. 1131 
The present NCSM/RGM resonance energies correspond 
to the energies where the diagonal phase shifts cross 90 
degrees. The calculation is as described in Table U and 
Fig. [5J We plot the absolute values of the energies as 



well as the excitation energies and the energies relative 
to the calculated and experimental thresholds. Overall, 
the NCSM calculation produces more binding by about 
600 keV as already discussed in the previous subsection. 
The NCSM/RGM generates excitation energies for the 
resonances systematically lower than the corresponding 
NCSM results. There is in particular a significant shift 
for the 2 + state. At the same time, as it can be seen 
from the bottom panel of Fig. [3j the excitation energies 
of the 2+0 and the l + states show a slower convergence 
rate with respect to the size of the HO basis expansion. 
This is a consequence of the inadequacy of the HO basis 
for the description of broader resonances. In this regard, 
the NCSM/RGM calculation is clearly superior. 



IV. CONCLUSIONS 

In this paper, we extended the ab initio NCSM/RGM 
approach to projectile-target binary-cluster states where 
the projectile is a deuteron. We gave details on the new 
aspects of the formalism and presented algebraic expres- 
sions for the integration kernels for the specific case in 
which the target wave functions are expanded in the SD 
HO basis. Among the new features, the dependence of 
the Hamiltonian kernel upon the three-body density of 
the target makes calculations technically challenging, due 
to the rapidly increasing number of matrix elements with 
the size of the multi-major shell basis. 

To test our formalism, we performed calculations for 
the bound and scattering states of the d-a system. In 
this case the three-body density calculation can be per- 
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FIG. 13: (Color online) Energies of 6 Li states calculated 
within NCSM and NCSM/RGM are compared to each other 
and to experimental data. Absolute energies (top panel), ex- 
citations energies (middle panel) and energies with respect to 
calculated and experimental d- 4 He threshold (bottom panel) 
are presented. Details of the calculations are given in Table [T] 



formed in a straightforward way using a closure relation. 
As the deuteron is weakly bound, its polarization and 
breakup cannot be neglected. A proper treatment would 
require the inclusion in the NCSM/RGM formalism of 
three-body final states: n-p-a. Although extensions of 
the approach in this direction are possible, this is quite 
challenging and has not been explored yet. In this first 
application, we approximated the three-body continuum 
by using deuteron pseudo-states. We compared our d- 
4 He results to experimental data as well as to a standard 
NCSM calculation for 6 Li using the same Hamiltonian. 
To facilitate benchmarking with the standard NCSM, we 
employed a soft SRG-N 3 LO NN potential with a low evo- 
lution parameter A = 1.5 fm -1 . In this way, we were able 
to reach convergence at N max w 12. The differences be- 
tween the NCSM/RGM and the NCSM results are then 
due to omitted correlations rather than to the adopted 
HO basis size. Interestingly, the NCSM calculation pro- 
duced a 6 Li g.s. energy lower by about 2%. This means 
that internal excitations of the 4 He target, neglected in 
our NCSM/RGM calculation, play some role and/or that 
the pseudo-state approximation of the three-body con- 
tinuum is not completely adequate. On the other hand, 
the NCSM/RGM calculation generates lower excitation 
energies for the broader resonances that present a slower 
convergence rate with respect to the HO basis expansion. 

Overall, the NCSM/RGM calculation is superior to 
the standard NCSM because it generates wave functions 
with proper boundary conditions for the bound state 
and, further, describes resonances and scattering states. 
However, to include all relevant excitations is a chal- 
lenge. Therefore, clearly the way forward is a unifica- 
tion of the two approaches. This can be accomplished 
by coupling the present NCSM/RGM basis, consisting 
of binary-cluster channels with just a few lowest excited 
states of projectile and target, with the NCSM eigen- 
states of the composite system as outlined in Ref. [3lj . 
Work on this unified approach is under way. 

Our immediate plans include the application of 
the NCSM/RGM formalism to the 3 H(c2,n) 4 He and 
3 He(<i,p) 4 He fusion reactions. This requires working in 
a NCSM/RGM model space including both ?i- 4 He (p- 
4 He) and d- 3 H (d- 3 He) channel states, that is a deuteron- 
nucleon (d, A^) transfer formalism which combines the 
deuteron-nucleus (presented here) and nucleon-nucleus 
(presented in Ref. 14]) formalisms as well as the inte- 
gration kernels resulting from the coupling between the 
(A— 1, 1) and (A— 2, 2) mass partitions, which will be the 
subject of a forthcoming publication. Our preliminary 
3 He(d,p) 4 He S-factor results were discussed in Ref. [32{ . 

The use of SRG-evolved NN interaction facilitates the 
convergence of the NCSM/RGM calculations with re- 
spect to the HO basis expansion. On the other hand, 
due to the softness of these interactions, radii of heav- 
ier nuclei become underestimated. Therefore, it is es- 
sential to further develop the NCSM/RGM formalism in 
order to handle thrcc-nuclcon interactions, both genuine 
and SRG-cvolution induced, in bound-state and scatter- 
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ing calculations. 

To apply the present deuteron-nucleus formalism to 
heavier target nuclei, i.e. heavy p-shell nuclei and be- 
yond, it becomes necessary to utilize the recently devel- 
oped importance-truncated NCSM dHHU. This gives us 
the ability to use large N ma _ x model spaces, that in the 
NCSM/RGM approach are of vital importance not just 
for the convergence of the target and projectile eigen- 
states but also for the convergence of the localized parts 
of the integration kernels [l5[ . 

Finally, our future plans also include a further gen- 
eralization of the formalism to projectile-target binary- 
cluster states with thrce-nuclcon ( 3 H, 3 He) and four- 
nucleon ( 4 He) projectiles. Calculations of the integra- 
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tions kernels for the three-nucleon projectile case are un- 
der way. 
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Appendix A: Hamiltonian kernel matrix elements 



Here we present matrix elements of the potential kernel (|22|) in the basis states (|18|) . For the first term on the rhs 
of Eq. J23]), we obtain 



U j : t 

SD\ k «*> 



Va-2,a-i(1 — Pa-2,a-i) 

1 



<I> 



SL> 



A - 



f h KI[ 1 f j b j' b K} f j' b j' c J bc \\T lT T{ 

2 hjtl bc Kr{ I' J I J I I'Ija J [jojb K 




T/,TT 7 J I T'T 2 1 



1 1 rp 

J-bc 



{a\,~a\ {KT) A-2axhT-\ y/1 + <W V 1 + ^ (b'c'J bc T bc \V\ bcJ bc T bc ) , (Al) 

SD \ V 1 I SD 



where we abbreviate a = n a l a j a \ etc. We note that the matrix elements of the interaction V in the antisymmetrized 

and normalized two-body basis are evaluated using just the first term of Eq. (j9|), i.e. Vij = VN(ij) + & ^^.J-^ 3 

(with Vn the nuclear part) as the average Coulomb interaction is taken care of with the help of Eq. (43) in Rcf 
For the second term on the rhs of Eq. (|23p we derive 



SD\ "*<■ 



Va-2,aPa-2,A-1 



SD 



hKvA r j'j h K\ r t iT t[ 




) V 1 + <W v 71 + S a ,b (a'c'lT 2 \V\abIT 2 ) . 

I SD 



(A2) 
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For the third term we get 



SD 



Va-3,a(1 — Pa-3,a)Pa-2,A- 
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E E E E 
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(A3) 



with the 12 — 7 symbols of the first kind [35| in the notation and the definition given in Appendix A of Rcf. [2C 
The matrix element of the fourth term on the rhs of Eq. (|23]l is obtained in the form 



SD 
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SD 
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Finally, for the last term on the rhs of Eq. (|23|) we find 
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